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Abstract. For lo-legged antiferromagnetic spin-1/2 Heisenberg ladders, a long-range spin pairing order can 
be identified which enables the separation of the space spanned by finite-range (covalent) valence-bond 
configurations into w + 1 subspaces. Since every subspace has an equivalent counter subspace connected 
by translational symmetry, twofold degeneracy, breaking translational symmetry is found except for the 
subspace where the ground state of «;=even belongs to. In terms of energy ordering, (non)degeneracy 
and the discontinuities introduced in the long-range spin pairing order by topological spin defects, the 
differences between even and odd ladders are explained in a general and systematic way. 

PACS. 71. 27. -(-a Strongly correlated electron systems - 75.10.Jm Quantized spin models 
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O 1 Introduction pointed out earlier in Refs. [^,@,|,|l for different sets of 

> 

■^ ■ long polymeric strips with graphite as the final member 

X; 

H , The discovery, about a decade ago, of high-Tc supercon- of these series, paralleling that of the square-lattice fam- 
ductivity ^ in hghtly doped "two dimensional" antifer- ily. Both theoretical and experimental studies ||,|l^,0|l|] 
romagnets and materials (initially) supposed to contain suggest that the nature of antiferromagnetic spin-i lad- 
coupled spin chains |@,|[|] , have generated a renewed in- ders with w=even legs differs from that of w=odd ladders, 
terest on low dimensional quantum spin-i systems. One For instance, w^even ladders are gapped systems, the gap 
of the concerns is the non-smooth crossover from one- vanishing exponentially with w, while u;=odd ladders dis- 
dimensional to two-dimensional systems (see, for instance, play characteristics similar to one-dimensional spin-i sys- 
Ref. [g and references therein). This fact has also been 
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terns, namely they are gapless, with a doubly degener- arbitraryly large but finite M. We refer to this space as 
ate ground state, breaking translational symmetry |10|,|l3|, model space, TV" . The reasonableness of Ti}" is based on 
n4]. Furthermore, spin defects are confined in ladders with the fact that the dimer-covering configurations (or Kekule 
w=even but they are not if w=odd. Numerical results |Q structures ||l^, as have always been termed in Resonance 
indicate that, in the infinite limit, the ground state of the Theory) are the lowest-lying monoconfigurational singlets, 
two dimensional system, towards u;=even and i(;=odd se- Thus, they provide a good zero-order picture. Then, on ap- 
ries must converge to, has long-range antiferromagnetic plying the Hamiltonian iJ, it can be noticed that the near- 
order and gapless excitations. ness of spin pairing tends to be preserved. When H is re- 
in this paper we will consider antiferromagnetic quan- stricted to the isotropic nearest-neighbour spin- 2 Heisen- 
tum spin-i ladders with w legs, (even) L ^ 00 rungs, berg Hamiltonian the shorter-range RVB picture should 
free boundary conditions in the inter-chain direction, and apply best for small even w, while i(;=odd or wide even w 
translational symmetry in the chain direction. It is as- ladders are expected to require long-range RVB pictures 
sumed that the Hamiltonian appropriate to describe these [QQ. For instance, M-range RVB pictures neglect cor- 
systems contains only short-range interactions preserving rections lying higher than the M order in Perturbation 
the total spin of the system. At half filling, we assume that Theory and have to be considered with caution. Never- 
the Hamiltonian which governs the lowest-lying region of theless, additional terms in the Hamiltonian, as frustra- 
the spectrum is the spin-^ antiferromagnetic Heisenberg tion, are expected to stabilise the finite-range RVB wave 
Hamiltonian, functions with respect to other Neel-based ansdtze (see 
rj _ v^ 7 c a (T\ [pO[ and references therein). In addition, there exist finite- 

range Heisenberg models for which short-range Kekule 



where S„i is the spin operator for spin on site ni, n in- 



dicating the rung and i the leg, and the Jni,mj are the 
exchange-coupling parameters. The Jni,mj £^re assumed to 
decrease very rapidly with distance, the nearest-neighbour 
Heisenberg Hamiltonian with isotropic J being the domi- In order to separate the model space TV" into non- 



structures are exact ground states and also short-range 
RVB ansdtze certainly apply for so-called "bond-dimer" 
models (see, for instance, [|l[|2|,|2|,||,||). 



nant part of H. mixing different subspaces, several attempts have been 

Since the ground state of such a Hamiltonian for a bi- made to find associated topological quantum numbers, 

partite system with equal number of sites in the two parts For instance, the occurrence of a topological long-range 

is known to be a singlet pGJ , resonating- valence-bond- type order (LRO) was first discussed p^ to rationalise the 

wave functions are defined in the space spanned by M- ground-state instability to bond alternation in spin-1/2 

range (covalent) valence-bond (VB) configurations, with linear Heisenberg chains. Simultaneously, this LRO has 
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also been discussed in the context of applications to con- The energy ordering among the lowest-lying state in 
jugated hydrocarbons [|6||7|,^,^. Latter, Klein et al. |p8| every subspace is estimate by the dimer-covering-coMntmf; 
and independently Thouless p9| introduced the gap or res- approximation P,B^p8] . Counting the dimer-covering con- 
onance parity, and Kivelson et al. [Q and Sutherland |Q, figurations has been achieved by a transfer-matrix tech- 
defined the winding number. These numbers allow the sep- nique (see, i.e., Refs. |34p5[ and references therein). Also, 
aration of the short-range VB states for odd-width strips since every subspace has an equivalent counter subspace 



in two subsets leading to degeneracy |p2,33|. The relation connected by translational symmetry, twofold degeneracy 

between topological LRO and winding numbers is given is naturally obtained except for the subspace including 

in Ref. |p3|. Also, in Ref. [ p3| , a resonance quantum num- the ground state of u'=even ladders, irrespective of the 

her, Dm which specifies the local (at boundary n) array details of any Hamiltonian preserving translational sym- 

of singlets, has been defined for VB systems with biparti- metry. Furthermore, in the present paper it is shown that 

tioning conditions. Still, arguments based on a topological a topological spin defect introduces a discontinuity in the 

LRO have been applied to the qualitative analysis of dis- LR-SPO, except for the ground state of u;=even ladders, 

tortious, excitations and their coupling for square-lattice Then, understanding energy ordering, degeneracy, and the 

strips p8| and, more quantitatively, to different polymers discontinuities introduced in the LR-SPO by topological 

p4,B5|. More recently, simple topological effects in short- spin defects, allow a general and systematic explanation 

range RVB were also predicated in Refs. |1^,Q for cou- of the differences between even and odd ladders. 

pled Heisenberg Chains, based on numerical results from 

This paper is organised as follows: In Sec. we show 
density matrix renormalization group (DMRG) techniques 

that the (covalent) VB configurations have a LR-SPO, 
on clusters. 

which allows the separation of the model space in differ- 
ent (asymptotically orthogonal and non-interacting) sub- 

Our purpose in this paper is to show that for anti- spaces. In Sec. the energy of the lowest lying state in 

ferromagnetic quantum spin-i ladders a long-range spin- every subspace is estimated within the dimer-covering- 

pairing order (LR-SPO) associated with the resonance counting approximation. In Sec. ^ we obtain the disconti- 

quantum number Z3„ can be defined. This LR-SPO allows nuity in the LR-SPO associated to the presence of a topo- 

to separate the model space Ti"^ into w + 1 subspaces. logical spin defect. In Sec. ||the results are presented and 

Configurations belonging to mutually different subspaces discussed. Finally, the conclusions are collected in Sec. ^ 
should differ repeatedly on each of the L rungs of the lad- 
der. Then, they are asymptotically orthogonal, and never 
mix by applying a few-particle operator. 
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2 Singlets and long-range spin-pairing order 2.1 Dimer-covering model-space approximation 

Quantum spin-i ladders with N = wxL sites, with (even) For any Kekule structure, let us define P+ (P") as the 

L -> oo and free boundary conditions along the interchain number of arrows pointing to the right (left) across a 

direction, are bipartite system with a singlet ground state, boundary, /„, lying midway between rungs n and n -f 1 

Therefore, the ground state can be written as a weighted (see Fig. |), and /„ as the number of SP with both sites 

superposition of a non-orthogonal complete basis set of ™ ^^e rung n, 

w — b 
singlets, | Si), i = 1 to dN, In — OA, ■ ■ ■ ^ , (3) 

TV! where 

"^""^ (iV/2-fl)!(iV/2)!- ^^^ fo,i« = even, 

b^l (4) 

In a bipartite system, sublattices A and B (starred) can [ 1: w; = odd. 

be identified and a set of cIn independent singlets can ^^ ^« (^n ) i^ the number of sites belonging to the inter- 
be constructed by pairing to a singlet each of the N/2 section of rung n and sublattice A (B), it can be easily 
spins in A to a spin in B. We represent one of these spin- seen tnat 

pairings (SP) by an arrow from the site in the sublattice A _ p- , p+ , t 

A to its partner in B (see, for instance, Fig. 1, where a r „ , „ 

'-' Wn ^ Pn-l + Pn + In. (5) 

complete set of linearly-independent singlets for six-site 

Choosing A and B sublattices according to 
systems are represented). Overlap, {st \ Sj), and matrix 

elements, (si I iJ I Sj) can be evaluated using the Pauling's Wp — Wq — b, (6) 

]l7|,p9| superposition rules. 
— — it can be written 

For the sake of simplicity, we first introduce the LR- 

u;;^ = i[u; +(-!)"&], 
SPO of VB configurations defining the (local, at boundary 2 



w, 



.B 1 



•^ = -[w-{~irb]. (7) 



n) resonance quantum numbers, I?„, when boundaries are " 2 

chosen to run parallel to rungs, and the model space is Substracting Eqs. (|) and using Eqs. (|) 

restricted to the dimer-covering approximation. Later we 

p-_,-p+_,+p+-p- = w^-w^ = {-irb. (8) 

show that the inclusion of longer pairings and/or using 

more general boundaries does not spoil this LR-SPO. The Defining the resonance quantum number D,, at boundary 

only effect of selecting boundaries of a different shape is •'" 

D = P"*" — P~ (9) 

changing the origin of the LR-SPO parameter. Finally, in " ~ " " ' ^ -' 

this section, we discuss the consequences of the LR-SPO we obtain 

on the eigenstates of H and their degeneracy. _D„ = Dn-i + (—!)"&. (10) 
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Then a SPO parameter D = Dq can be associated to any states allow sites in v4-sublattice to be SP to sites in B- 



dimer-covering configuration, so that 



Dn = D-^[l-{-ir]b. 



(11) 



Since 



D can take w + 
D = 



P+ = 1 w^ 

p,7 = o,i,...w^, 

1 different values, 
w + b w + b — 2 



(12) 



b- 



(13) 



sublattice no more than 3 bonds apart. These states can 
be directly generated by the "re-coupling" Q,^ of two 
simply neighbouring dimers, i.e. unlinked bond-pairs with 
one and only one site in a pair being a nearest neighbour 
to a site in the other pair. It is worth noting that these re- 
couplings satisfy the non-crossing rules. Then, the 3BR-SP 
model subspace incorporates any singlet obtained from a 
dimer-covering singlet allowing an arbitrary number of un- 
linked re-couplings of two simply neighbouring bond-pairs. 
Still longer-range model spaces can be obtained allowing 
and the (dimer-covering) model space can be partitioned 5BR-SP, 7BR-SP, . . ., to be included in W'. Nevertheless, 
in w; + 1 subspaces, H^, according to the value of D. singlets with very long bond-range SP should contribute 

less, so a reasonable model space will be that including 
2.2 General formulation of the LR-SPO singlets with SP up to M bonds apart, M not necessarily 

small. 
The LR-SPO and D introduced above have been related, 

r,, , r-i--, i-iij- • Let us now also allow the boundary f^ to be a line 

tor the sake of simplicity, exclusively to dimer-covermg ^ Jn 

n .. J J- u J • ■ 11 1 i running from one side of the ladder to the other side, with 

conngurations and to boundaries running parallel to rungs. '^ ' 

^T 4-u i • J-- 117 -11 j-u J- j-u nl being the first site to the left of f,f in leg 1. We assume 

Now we remove these restrictions. We will see that the ^ Jn t> 

, r ii u J • 1- -i- r J- r ii 1 J 1 that f,? Can go up and down, but it is self-avoiding and is 

shape of the boundaries limiting fragments of the ladder Jn b t- ; & 

,, . 1 . r 1 ■ • • ■ ■ 1 i not hitting any site. Thus, f^ must unambiguously break 

or the inclusion of long-range spm-pairmgs is irrelevant o j ^ ju o j 

and that the LR-SPO can stiU be defined. ^P *^^ ^^^^^"^ ''' *^° '"^S^"^'^ ^^^* '"^gi"^' ^"' ^^^ "S^* 

mi ,. • 1 , . ^ . ■ , 1 region, i?„ (see Fig. H). Therefore, two non-intersecting 

i he dimer-covermg model space is not invariant under o ' " \ o Ijy ' o 

4-1, XT -14- • 4- T7 • 4- 4-1, w 4- boundaries, /.f and /^, , n < m, define a fragment, F^„, , 

the Hamiltoman operator, i^or instance, the XY terms, " ™ ^ "'™ 



'^ni'^mi' of the nearest-neighbour Heisenberg Hamiltonian 



as the intersection of -R„ and L,, 



acting on a Kekule structure yield singlets with SP be- We define P,f+ (-Pn ) ^^ the number of arrows pen- 

tween sites up to 3 bonds apart (see Fig. |3|). Then, as etrating the boundary /,^ with the arrowhead in the -R„ 

a first step, linearly independent singlets with 3-bond- (Ln) region. /,^ „j is the number of arrows with both ends 

range (3BR) SP should be incorporated into H^ to go in F,f „. ^^ „ (''nm) i^ the number of arrows starting in 

beyond the dimer-covering approximation. These 3BR-SP -/?,„ (Ln) and with the arrowhead in -L„ (-/?,„), i.e. with 
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no partner belonging to F^ „. Finally, F^"^ (^n m) i^ the and f2j do not intersect. Ip is the number of SP with both 
number of sites belonging to the intersection of F-( ^^ and sites in Fq2j', If ('"f) is the number of arrows connecting 
sublattice A (B). Then an A site in i?2j (^o) to a S site in L^ (i?2j), and Ff^2i 



^9 ,7-5 ^ f^gA ^" o.2j 

''n.n-\-'p ' ^ n+p n,n+p ^ n.n-\-p n,n-\-p^ 

n.7i+p ' n+p n.7i+p ' n.7i+p n,n+p'\ J 



Subtracting these two equations we obtain 



^n+p ^n — ^n,n+p ^ n,n+pJ V-^^7 



{F^2j) is the number of sites belonging to the intersection 
of -fo,2j and sublattice A (B). Then 

P(?- -If + P+ -rF+lF = F„% , 



Pi^ -rF + P^j -If+If^ <2, • (20) 

Subtracting these two equations 



with 

Di = />,?+ - P9- ^ (16) - D^ + D2j = F^2j - Fi2j ■ (21) 

Let us analyze -F^_„+p leg by leg. From Fig. |, it is Using Eq. (|l|) we obtain 
readily seen that 



Da = D- {F,% Fi,^ , (22) 






which is independent of j, provided it is not too small to 
g 



^n,n+2j+l ^n,n+2j + l — %' \^ ' ) 

prevent /q and f2j are intersecting. Thence, there is a one 



where j is a positive integer with the restriction /» and to one correspondence between allowed values of D^ and 
fn+2j "io i^ol intersect, and &f^ = fo^_|_2 for any n. Then, £) 
choosing n — 0, 

n9 n9 — n ^•^ LR-SPO of eigenstates and degeneracy 

^2j ^ ^0 — U, 



^2j + l ^0 — "O- 



^ I So far, we have separated the model space in w + 1 sub- 



, • .; r 1, ,, , m-,rA . 7->n 7->o i spaccs. At thls point we note that two singlets from dif- 

Agam it follows that a SPO parameter D^ = L)g can be 

. ^ 1 , ,rn n i- ferent subspaces must be different repeatedly at every po- 

associated to any VB conhguration, so ^ i- j j f 

sition along the ladder. Therefore, according to Pauling's 

^^ = ■0^- J[l-(-l)"]&^ (19) ^r^ 

^ island-counting technique [g/7|,|39|] , they are asymptotically 

with b9 = bg. orthogonal and non interacting via any interaction medi- 

The general order parameter D^ can be related to the ated by a few-particle operator. Then the matrix of the 

previous one, D, in a simple way. For simplicity, without Hamiltonian asymptotically block-diagonalises, so config- 

loss of generality, let us consider the fragment ^0 2^ liniited urations belonging to different subspaces do not mix in 

by /q and f2j, j > 0. We select j in such a way that as /g the configuration-interaction sense. Thus D may be taken 
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as a long-range order parameter labelling the eigenstates for a class of benzenoid hydrocarbons M (with C=0.5667) 

^D of the D block. and for finite square-lattice fragments |Q (with C~0.75), 

Furthermore, 7i^ and 'H^_q subspaces are equivalent by fitting the logarithm of the Kekule-structure count to 

by a translation T of one rung along the legs. Then, the resonance energy calculated from an equally- weighted 



T^D = ^b-D (23) 

Consequently, degeneracy is always expected to occur, ex- 
cept for D = Q and w=even, i.e. 6 = 0. 



Kekule-structure wave function. 

Arguments supporting Eq. (|2J) relay on the fact that 
the energy is an extensive magnitude, i.e. scales as the sys- 
tem size, while the number of Kekulc structures needs to 



scale exponentially with the system size. Then, C may 
6 energy ordering l^g determined by fitting the logarithm of the Kekule- 



Within the dimer-covering approximation the resonance 
energy in units of J, Ner{w, D), i.e. the ground-state en- 
ergy correction below the energy of a single Kekule struc- 
ture (-0.375-/V), depends on the configuration interaction 
amongst the different Kekule structures. It has been ar- 
gued P,p8[ that one might consider this energy lowering 
to depend solely on the dimension of the space spanned 
by the appropriate Kekule structures. Let n{w, D) be the 
dimension of the D block for a w-legged ladder. Since 
n(w, D) is multiplicative in terms of a break up into sub- 
systems while the energy is additive, such a functional The values of n{w, D) can be easily obtained (see the 



structure count to better estimates to the energy than 
those used in Refs. [§ and |3^. As a first attempt to 
do so, we obtain C^=0.84, 1.0 and 0.93 appropriate for 
square-lattice strips with L —i- oo and w=2, 3 and 4, re- 
spectively, when the "exact" energy values of Table II in 
Ref. ^^ have been used. The weighted average of them 
could be used as a rough estimate of C appropriate for the 
two-dimensional square lattice. Assuming that the error of 
Cw typically behaves as 1/w, we obtain C = 0.94 ± 0.19, 
although a more reliable value would be desirable. 



dependence should be of the form Appendix) as the L/2 power of the highest eigenvalue, 

„ J ^wD' ^^ t^^ -^ block of the square of a dimer-covering- 

er{w,D)^ -\nn{w,D), (24) 

^^ counting transfer matrix T^ . Then, the resonance energy 

where C is a fitting parameter independent of the struc- of the lowest lying state in the subspace D can be written 

ture to some degree, and in particular presumed to be (at as 

least asymptotically for large w) the same for all w-legged 

CJ 

ladders. For energy-ordering purposes the value of C is ir- ^ w " ' 

relevant. Nevertheless, for the nearest-neighbour isotropic 
Heisenberg model the value of C have been determined 
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4 Topological spin defects and LR-SPO Furthermore, it is possible to form a local region be- 

tween rungs n and m of LR-SPO D ± 1 placing two spin 

Let us now suppose that there is a topological spin defect 

defects in a wave function of LR-SPO D, one in a site ni 

on site i of rung n, i.e. for any reason the site ni remains 

and the other in a site mj, with n + i + m + j=odd. 

not spin-paired, 



ni £ { 



A, n + i = odd, 

(26) 5 Results and discussion 

B, n + i = even. 



Then, choosing boundaries parallel to rungs, Eqs. (0) be- Considering that u;-legged (even) i -> oo antiferromag- 

netic spin-i Heisenberg ladders are bipartite systems with 

singlet ground state, we have assumed that their ground 

^ state is described by an RVB many-body wave-function 

" " ^ " "■ 2 ^ ^ ansatze withm the space spanned by the M -range (co- 

and the recurrence relation across the rung n with a topo- ^^l^'^*) valence-bond configurations, H^ , with arbitrarily 

logical spin defect on site ni is rg ^ . 

We have shown that any (covalent) VB configuration 

A. = A.-1 + (-1)" [h + (-1)'] . (28) 

exhibits a LR-SPO. This LR-SPO is characterised by a 

If the order parameter to the left of the spin defect is Di, parameter D, which can take w + l different values, and 



according to Eq. ([TT]), the order parameter to the right of allows to specify the local (at boundary) array of spin- 
site ni, Dr, will be pairings penetrating the boundary n, for any n. The shape 

of the boundary selected to define D is quite arbitrary: it 
Dr = Di + [-lT+\ (29) 

can go up and down, but it is self-avoiding and is not 

Then, a topological spin defect can be seen as a domain hitting any site, so it is able to break up the ladder in two 

wall separating sectors with order parameters Di and Dr. (non- longitudinal) parts. When the boundary is chosen to 

Furthermore, since the subspaces D and 6 - Z? are de- be parallel to the rungs, the allowed values of D are: 

generate, the energy per site associated to the sectors to nn^h in^h~9 h — 

the right and to the left of a spin defect located at the ni 

site will be degenerate if where b is zero (one) for it;=even (odd) . Choosing another 

boundary, the values of D simply change by a boundary- 

D, = -[b^i^ir+^] (30) 

^ dependent integer. 

This equation has a solution only when 6 = 1, i.e. w is This long-range spin-pairing-order allows to separate 

odd, and Di = (with D,. = 1) or Di = 1 (with Dr = 0). TC" in w+l subspaces, TCfj- Two configurations in different 
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Hf) differ repeatedly on each of the L rungs of the ladder, energy ordering 

Consequently, they are asymptotically orthogonal and non 

Eb<Eb+i<---<Ei^^^,y (34) 

interacting via any interaction mediated by a few-particle 

operator. Thence, the matrix of the Hamiltonian asymp- Therefore, for 'u;=even, the ground state belongs to the 

totically block-diagonalises. Thus, the eigenstates oi H do non-degenerate D — subspace, with 

not mix configurations with different value of D. There- 

T1'o = ^o, (35) 

fore, D may be taken to be a long-range order parameter 

labelling the eigenstates of the D block. This gives rise to while, for i(;=odd, the ground-state manifold is spanned 

w + 1 RVB Anzdtze exhibiting different LR-SPO, 'Pd- It by fFo and tf'i, which are eigenstates of H, but not of T. 

is also worth noting that asymptotic strong orthogonality The eigenstates 'P± of the translation operator, are defined 

and LR-SPO in RVB wave functions persist with a long in the ground-state manifold, 

finite cut-off. -i 

<P± = -= (tf'o ± ^i) , (36) 

Furthermore, subspaces 7i^ and 'H^'_jj are equivalent 

by a translation of one rung along the legs, T. Then so 

T<P±=±<P±. (37) 

T^D=^b-D, (32) 

with fc == and k — t:. This RVB asymptotic degeneracy 

and the energy of the corresponding wave functions must for w^odd is consistent with a very wide body of evidence. 

be It could be argued that the true ground state of the 

spin-i nearest-neighbour Hcisenberg Hamiltonian for w = 
Ed = Eh_D. (33) 

1 and L — > oo should be described in terms of two degen- 

Thus degeneracy is always predicted to occur, except for erate subspaces, in order to regain the non-degenerate and 

u;=even and D = 0, irrespective of the details of any gapless Bethe Ansatz state. For finite (even) L with cychc 

Hamiltonian preserving translational symmetry boundary conditions the two degenerate RVB wave func- 

We have obtained the values of A^^iu of Sec. for tions do interact, in the configuration interaction sense, 

w ranging from one to twelve and D from 6 to {w + giving a splitting of the translationally-adapted states ^± 

6)/2. Table H summarises the zero order resonance ener- (with k — and k — tt), which leads to a non-degenerate 

gies Eriw^D), in units of C, calculated by using Eq. (P5|). ground state. Thus, for finite L, restriction to just one 

Since it is unlikely that the zero-order energies are dras- subspace does not yield a true ground state. Nevertheless, 

tically modified by the small corrections to the isotropic when the splitting is examined as a function of L, it is 

nearest-neighbour Hcisenberg Hamiltonian, we expect an seen that the gap closes at least as fast as 1/L and the 
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states ^± approach to an accidental degeneracy. This fea- two sectors of the ladder with order parameters Di and 

ture is basically related with the dimerzing instability of Dr 

the Bethe ansatz state. For instance, if a small dimerizing jj^ = Di + (—1)"^' (39) 



Furthermore, the energy per site associated to the right 
and to the left sectors of the spin defect will be degenerate 
only when 6=1, i.e. w is odd, and Di = (with D^ = 1) 
or Di ^ 1 (with Dr = 0). Also, it is possible to form a local 
region between rungs n and m of LR-SPO D ±1 placing 
two spin defects in a wave function of LR-SPO D, one in a 
site ni and the other in a site mj, with n + i + jn + j=odd. 



interaction occurs, then the accidental degeneracy of Eq. 
( p3| ) no longer holds. Therefore, either tf'o or ^i will be 
lower in energy and completely dominate the wave func- 
tion, leading to a spin-Peierls broken-symmetry ground 
state ||. 

The w — > oo limit is obtained by fitting Sr{w, 0)/C by 
a power series in 1/w, 

£r(-W,0) _ £r(oO,0) , ^ , _^ /3g^ 

C ^ C w w"^' Now, half-filled excited states (other than the d„ — 1 

InTablegthe values of £r(oOi 0)/C, fli and a2 are collected singlets orthogonal to the ground state) or even slightly 

for both w=even and w=odd, along with their weighted doped states are analysed via topological spin defects, 

average for £^(00, 0)/C'. For the purpose of seeing how rea- There are different types of excitations conceivable from 

sonable this simple model is, a value of C is needed. When a Maximally-spin-paired ground state. For instance, pre- 

C is fitted to dimer-covering estimates of the energy for serving half filling (one electron per site), there are pri- 

square-lattice fragments, the value 0.75 is obtained ||38| . marily spin excitations. In this case, two topological spin 

Using this value of C, the present simple model predicts an defects, one in an A site and the other in a _B site, are ob- 

energy per site of ~ —0.594 (in unit of J), in good agree- tained by breaking one SP to form a triplet state. Away 

ment with the value —0.604 calculated at the same level from half filling, removing (adding) one electron produces 

by Liang, Doucot and Anderson Efl] . When C is evaluated a vacant (doubly-occupied) site and the ladder becomes 

as the weighted average of the Cw, w=2, 3 and 4, derived a doublet, with two sites, again one in an A site and 

from the "exact" estimates of the energy of Zivkovic ef aZ. the other in a B site, that cannot be SP. In this case 

|p8| (see section^ we obtain C = 0.94±0.19. Obviously, a hopping terms must be retained in the Hamiltonian and 

more reliable C value would be desirable. Using this value the so-called t-J model applies. Thence, the doublet is a 

of C, an energy per site of w —0.65 ± 0.10 (in unit of J) weighted superposition of configurations with a single spin 

is obtained, which also compares fairly well with the best not SP and a vacant (doubly-occupied) site. Either a va- 

estimatc (—0.668) of Liang, Doucot and Anderson E^]. cancy or a doubly-occupied site may also be assimilated 

We have shown that a topological spin defect, i.e. a to a topological spin defect, although there is no spin asso- 

non-paired site, can be seen as a domain wall separating ciated with them. Then, away from half filling, it may be 
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conceivable a local region limited by two vacant (doubly- pendently. Consequently, a gapless triplet spectrum is not 
occupied) sites, or even a vacant (doubly-occupied) site inconsistent with the results of this work. This feature can 
and a single non-SP spin (provided that the doping is not be understood as a generalisation to any odd-legged spin- 
so strong as to preclude a maximally-spin-paired ground 2 antiferromagnetic ladders of Lieb, Shultz, and Mattis 



state). Therefore, there are low-energy spin and charge ex- theorem [|13, 14| holding for one-dimensional systems. Sec- 
citations. Still, going up in the hierarchy of Hamiltonians, ond, away from half filling, removing (adding) an electron 
the Hubbard or even a more general Hamiltonian has to to the system yields a non-confined pair of sites, one being 
be considered. In this case, still another type of excitations a vacant (doubly-occupied) site and the other a non-SP 
(though presumably of higher energy if a Heisenberg-like site. It is worth noting that the vacant (doubly-occupied) 
Hamiltonian is assumed to govern the lowest-lying region site holds the charge, while the non-SP site holds the up 
of the spectrum) can be produced relaxing the single- or down spin, leading to charge-spin separation. 

occupancy constrain. This leads to the ionic states, i.e. 

In clear contrast, when w=even, the order parame- 
states with at least a couple of sites, one doubly occupied 

ter of the local region limited by the couple of spin de- 
and the other empty. Therefore, the couple of spin defects, 

fects is always associated to higher energy per site. This 
associated to an excitation above a Maximally spin-paired 

indicates that the couple of spin defects should remain 
state 'Fd, are limiting a local region of LR-SPO D ±1. 

as close as possible so confinement is predicted to oc- 

When w=odd and £'=0 {D=l), there can be local re- cur. Hence, at half filling, the energy difference between 
gions with D = ±1 (-0=0, 2). Then, since subspaces D—0 the triplet and the ground state will be finite, and the 
and D—l are degenerate, it is possible to have a local re- w=even ladders are expected to be gapped. Neverthe- 
gion with identical per site energy inside and outside the less, we observe that the energy difference per site be- 
local region. In this case the topological spin defects lim- tween the lowest lying subspaces, Ae, behaves as ~ w~", 
iting the local region are not confined, though it may hap- with a sa 1.8, i.e. Ae decreases faster than 1/w. There- 
pen that they attract one another (with an ordinary short- fore, a lowering of confinement and the closing of the gap 
range potential). Two conclusions can be drawn from this is predicted for increasing values of w. At this point, it 
result. First, for the half filling case, triplets with the two could be argued that an excitation energy lowering can be 
spin defects very far apart from one another are possible, achieved by allowing a Bloch superposition of analogous 
Although breaking a singlet does cost some energy (due to couples of topological spin defects. However this lowering 
contribution of the diagonal terms), there is a gain in ki- is not expected to be as important as to close the gap, at 
netic energy (off-diagonal terms contributions), since (for least for w—2 [pSi. Nevertheles, within the scope of the 
i(;=odd) the two spins arc not confined and can move inde- present work it is not able to predict the energy order- 
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ing of the lowest singlet and triplet excited states. Now, bond, while a short-range RVB is expected for w=even 

away from half filling, when removing (adding) one elec- ladders ||l8|,[19[. 

tron the vacant (doubly-occupied) site will be bound to The above results also apply to open boundary con- 
the non-SP spin, forming a "quasi-particle" with charge ditions along the L direction. In this case, no pairing is 
and spin. Again we expect the lowering of the confinement expected to the left of rung n = 1, and to the right of 
and the charge-spin separation with increasing values of rung n = L. Then, open boundary conditions determine 
w. This quasi-particle is far from being a static vacant the long-range order parameter to be I? = 0, with a pair 
(doubly-occupied) site bound to a static non-SP site. The of spin defects still delimiting a, D — \ fragment, which 
hopping term of the t-J Hamiltonian allows the vacant is degenerate with that of 13 = solely for w=odd, so 
(doubly-occupied) site to move while the exchange part results above still apply. Furthermore, if the boundary 
mixes up all the doublet VB configurations with the con- conditions are so that would force an order parameter of 
strain of keeping the vacant (doubly-occupied) site and energy £r{w,D) {er(w,b~ D)), D ^ b, we expect that a 
the non-SP spin close to one another. Therefore, the ap- state with D — b+1 spin defects located close to each end 
propriate wave function must be a weighted superposition of the ladder, these defects limiting an extended area of 
of all possible static configurations in order to benefit from resonance energy per site er(w, b), would be favoured, 
the hopping contribution to the energy, while lowering the The results of the present work relay on two conditions, 
exchange energy by keeping both the spin and the vacant First, the bipartitioning condition, which allows to estab- 
(doubly-occupied) site bound. When removing (adding) lish the LR-SPO. Second, the energy ordering given by 
two electrons, arguments based on the LR-SPO alone are Eq. B3. This energy ordering is valid when the dominant 
not able to decide if the two-holes (two-electrons) state part of the Hamiltonian is the isotropic sin-1/2 Heisenberg 
is described as two quasi-particles or two bound vacant Hamiltonian. It might happen that the bipartitioning con- 
(doubly-occupied) sites. Nevertheless, we expect the two dition is not fulfilled and/or the corrections away from the 
vacant (doubly-occupied) sites to be confined to benefit isotropic spin-1/2 Heisenberg Hamiltonian are as impor- 
from the energy lowering due to the hoping term of the tant as the energy ordering is not that of Eq. (pj) , then 
Hamiltonian, as has been pointed out in Ref. |3g] , using the above discussion will not apply, 
numerical results from DMRG techniques on clusters. An- 
other argument to take into account is the range of the 5 Conclusions 

RVB. Since the two quasi-particles are expected to couple 

Identifying a LR-SPO and understanding energy ordering 
to a singlet, non-bound quasi-particle would imply a long 

and degeneracy among RVB wave functions, as well as the 

discontinuities in the LR-SPO introduced by topological 
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spin defects, allow a general and systematic explanation of arc 2™ different local states for each boundary, which can 

the differences between even and odd ladders. It provides be classified according to the value of D„, | e^j). The 

an intuitive understanding of the physics of spin ladders local states of position n -I- 1 are mirror images of those of 

in general. In particular, for odd leg ladders position n. 



1. a doubly degenerate ground state, breaking transla- 



A diTaer-covermg- counting matrix, Tn, is defined as 



Then, the number of dimer-covcring states in a 7i^ sub- 
space is 



tional symmetry, is obtained; I 1, | e„,/) can succeed | Cn-ii) 

(e„-i7 I T„ I e„j) = < 

2. a gapless triplet spectrum is consistent with the results I 0, otherwise 

(40) 
of the present work; 

3. a simple understanding of the deconfinement of spin 
excitations as well as charge-spin separation is pro- 
vided. n{w, D) = Y^ie^i I ^1^2 ■ ■ -T^ \ e^j). (41) 

On the contrary, for small w=even ladders n- r i- • • 1 t-, rr rr 

Since for any dimcr-covering singlet Dn-i = Dn+i, Tn%i+i 

1. a non degenerate ground state is obtained; is a block-diagonal symmetric matrix that does not de- 

2. triplet spin excitations are found to be gaped; pend on n (apart from the direction of the arrows in the 

3. a simple understanding of the confinement of spin ex- local states that it relates) we can omit the sub- index. For 
citations is provided; L ^ 00, the highest eigenvalue A^j-, of the D block T^ 

4. a lowering of the confinement and the closing of the dominates, and 

gap is predicted for increasing values of w. n{w, D) « A^d- (42) 

The author acknowledges valuable discussions with D. J. Klein. 
This research was supported by the DGICYT (project PB95- 
0884). 



Dimer covering counting 

Let us analyse from a local point of view the dimer-cover- 
ing singlets. We can identify a dimer-covering local state, 
I e„/), / ranging, according to which legs have an arrow 
across the /„ boundary. The direction of any arrow is fixed 
by n and the leg number. Then, it can be seen that there 
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Figure captions 

Fig. 1. A complete set of linearly independent (covalent) sin- 
glets for a six-sites system: a) 1 x 6; b) 2 x 3. 

Fig. 2. A fragment of a w—9 ladder, showing (inside a cir- 
cle) the value of D„ at each boundary /„. Note that D„ — 
Dn-i — ( — 1)"6 (sublattices A and B are selected as w„ —w^ ~ 
(-l)"b). 

Fig. 3. A dimer-covering configuration is an eigenstate of 
Sj • Sfe when sites j and k are spin-paired. The "ofi'-diagonal" 
singlet with {i,l*) and {k,j*) pairings is also obtained when 
Sj ■ Sfe acts on a VB configuration with spin-pairings {i,j*) 
and {k, l*). 



Fig. 4. A fragment of a w=7 ladder, showing (inside a circle) 
the value of D„ at each boundary /„. Note that still D„ — 
Dn-i + ( — 1)"& (sublattices A and B are selected as w„ —w^ — 
( — 1)"6), although here there are longer-range spin-pairings. 

Fig. 5. Portion of a ladder showing boundaries /^ and /^ 
which go up and down: a) m — n=even; h) m~ n=odd. These 
non-intersecting boundaries define a fragment of the ladder 
(the set of sites with a full dot or star) , f ^.^ i as the intersection 
of Rn and Lm ■ Analysing the fragment a leg at time, it is worth 
noting a series of facts: First, there are m — n sites in each 
leg. Also, whenever a site qi belongs to -F^_m, neither the site 
T"'~"qi nor the site T^'^^'^qi {qi translated m-n steps to the 
right or to the left) are included in the fragment. Furthermore, 
in (a), T"^^"qi remains in the same sublattice than qi while, on 
the contrary, in (b), T^^^^qi belongs to a different sublattice 
than qi. Finally and more important, in (a) half of the m — n 
sites belong to the sublattice A and the other belong to the B 
sublattice, while in (b) the fragment has different number of 
sites in the sublattice A than in the sublattice B. 

Table captions 

Table 1. The absolute value of the resonance energy in units of 
C, (In Ad)/'w, for the lowest lying state of subspaces with order 
parameters from D — to (w + b)/2, and the extrapolations 
to «; — ^ oo for the lowest lying subspace of both, the w=even 
and ui^odd, series. 
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Table 2. Parameters of the Eq. (Bq), which fits the ground- 
state (D — 0) resonance energy, in units of C, for the u;=even 
and ui^odd series and their weighted average. 
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